Effect of the Wyckoff position of K atom on the crystal structure and electronic 

properties of KFe 2 Se 2 compound 
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By the first-principle electronic structure calculations, we study the effect of Wyckoff position 
of K atom on the crystal and electronic structures of KFe2Se2 compound. When the K atoms 
take up the Wyckoff position 2a, 26 or 4c (the related structure of KFe2Se2 is called as Struc-2a, 
Struc-2b or Struc-4c), the calculated lattice constant c is in the range of 13.5~14.5 A, 15.5~16.7 
A or 18.6~19.1 A. Three concentric cylinder-like Fermi surfaces emerge around T-Z in Brillouin 
Zone for the Struc-2b in nonmagnetic state, different from the Struc-2a and Struc-4c. The Fe-Se-Fe 
angle is 107.8°, 108.8° or 110.7° in the collinear anti-ferromagnetic state and the superexchange 
interaction J2 between two next neighbor Fe moments is 13.08 meV/S 2 , 20.75 meV/S 2 or 11.86 
meV/S 2 in Struct-2a, Struc-2b or Struc-4c structure respectively. The Struc-2b and Struc-4c have 
a good correspondence to the newly discovered superconducting phases with T c =40 K and T c =30 
K in KFe2Se2. Our findings provide a reasonable approach to understand the existence of multiple 
superconducting phases in alkali metal intercalated FeSe superconductor. 

PACS numbers: 74.70.Xa, 74.20. Pq, 74.20.Mn 



INTRODUCTION 

The recent observation of superconductivity about 30 
K in the layered iron selenides compound Ko.sFe2Se2 [lj 
initiated extensive research on the iron selenide supercon- 
ductors, which was formed by intercalating alkali metal 
between FeSe layers and had the ThCr 2 Si2-type struc- 
ture (Fig. 1(a)). In contrast to the iron arsenide su- 
perconductors, A !/ Fe2- a; Se2 (A = K, Rb, Cs or Tl) is 
a new class of iron based superconductor. Their elec- 
tronic structures are very different from those in the iron 
arsenide compounds, that is, only electron Fermi sur- 
faces are presented around the zone corners while no hole 
Fermi surface near the zone center, which is confirmed by 
the experimental measurement and the previous first 
principle calculations Q. 

The crystal structures of Aj,Fe 2 _a;Se 2 are complex 
since for different x values there are different Fe vacancy 
arrangement schemes. The X-rays, transmission electron 
microscopy, and neutron scattering measurements indi- 
cate that the composition of K-intercalated FeSe super- 
conductors was close to K !/ Fei.6Se2 with a fivefold ex- 
pansion of the parent ThCr 2 Si2 unit cell in the a&plane, 
namely a \/E Fe- vacancies superstructure 0,0]. An- 
other experiment show that there exists a rhombus Fe 
vacancy pattern related to K y Fei.5Se2 @. However it is 
reported that in the same sample there exist two distinct 
phases, the insulating phase with well-defined y/S x Vb 
Fe- vacancies superstructure and another KFe2Se2 phase 
containing no Fe vacancies, the latter being suggested to 
be responsible for superconductivity The composi- 
tion and structure responsible for the suerconductivity is 
still in debate. 

Most recently, Li, Na, Ca, Ba, Yb, or Eu layer is in- 



serted between adjacent FeSe layers to optimize the su- 
erconduting properties of the iron selenides compounds 
It is the surprising phenomenon that two supercon- 
ducting phases were observed in Ko.sFe2Se 2 with the one 
T c of 30 K and another T c of 40 K. In addition, two 
kinds of crystal unit cell with unexpectedly long lattice 
constant c=16.14 A and c=20.48 A were found to corre- 
spond to two superconducting phases, and the c=20.48 
A unit cell disappeared after one hour's exposure to air 
accompanying the T c change from 40 K to 30 K. In the 
other reports, two superconducting phases were ever ob- 
served in Ka;Fe2-ySe2 [§-11 • However, there are rare 
theoretical studies on the difference of crystal and elec- 
tronic structures between these superconducting phases, 
and the reason for higher Tc in new superconducting 
phase. 

The Wyckoff position of K atoms has a great influence 
on c axis length of crystal unit cell and the Fermi surface 
characters in KFe2Se2- In experiment the occupations of 
different Wyckoff positions by alkali metal atoms are ob- 
served [HI]. Meanwhile, the fact that no clear peak due to 
N-H vibrations is detected in infrared spectroscopy mea- 
surements |8[ excludes the existence of NH3 in KFe 2 Se 2 
synthesized by the liquid ammonia method. So we can 
deduce that three superconducting phases of KFe 2 Se 2 
with different lattice parameter c should correspond to 
the different Wyckoff position occupied by K atoms. In 
the paper, we focus on the effect of different Wyckoff po- 
sition of K atom in space group I4/mmm on the crystal 
structure and electronic properties of KFe2Se 2 , with K 
atom sited in 2a, 2b or 4c Wyckoff position respectively. 
The paper was organized as follows. Firstly, we optimize 
the crystal structure and obtain the reasonable lattice 
constants of KFe 2 Se 2 , in agreement with the experimen- 
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FIG. 1: (Color online) Schematic structure of the tetragonal 
unit cell of Ka;Fe2Se2 containing two formula units: (a)K atom 
is situated in 2a Wyckoff position, (b)K atom is situated in 26 
Wyckoff position, (c)K atom is situated in 4c Wyckoff position 
and only two of four 4c positions are occupied, the a and Hse 
indicate the Fe-Se-Fe bond angle and the Se height to Fe 
plane. 



tal values of c=14.04 A, c=16.16 A, and c=20.48 A in 
three superconducting phases. Secondly, the electronic 
structures of Struc-4c is similar to those of Struc-2a im- 
plying they have the similar superconduting behavior. 
The Fermi surface of Struc-2b exhibits the different char- 
acters, that is, the concentric cylinder like Fermi surfaces 
are presented around T-Z, which should have a neces- 
sary connection to the superconductivity with the higher 
Tc=40 K. Lastly, from three aspects of the Fe-Se-Fe an- 
gle, electronic structures and the next neighbor superex- 
change interaction J2, we demonstrate that the Struc-2b 
structure is closely related to the superconducting phase 
with T c =40 K, while Struc-4c to the superconducting 
phase with T c =30 K in KFe2Se2 synthesized by the am- 
monothermal method. 



METHOD AND DETAILS 

In our calculations the plane wave basis method was 
used [13|. We adopted the generalized gradient approx- 
imation (GGA) with Perdew-Burke-Ernzerhof formula 
14j for the exchange-correlation potentials. The ultrasoft 
pseudopotentials were used to model the electron-ion 
interactions. After the full convergency test, the kinetic 
energy cut-off and the charge density cut-off of the plane 
wave basis were chosen to be 800 eV and 6400 eV, re- 
spectively. The Gaussian broadening technique was used 
and a mesh of 16 x 16 x 12 k-points were sampled for 
the Brillouin-zone integration. In the calculations, the 
lattice parameters with the internal atomic coordinates 
were optimized by the energy minimization. The Struc- 
4c tetragonal unit cell should have included four K atoms, 
but only two K atoms is included in order to keep the 
same composition K:Fe:Se=l:2:2 for Struc-2a, Struc-2b 
and Struc-4c. 



RESULTS AND DISCUSSION 

For the first part, we study the structural proper- 
ties of KFe2Se 2 including the lattice parameter, the Se 
height to Fe plane Hg e and the Fe-Se-Fe bond angle a in 
FeSe layer. Alkali metal intercalated iron selenides have 
the tetragonal crystal structure with the space group 
I4/mmm. Fe atom and Se atom are sited in the 4d and 
4e Wyckoff position in this space group respectively, con- 
structing the stable FeSe layer in KFe2Se2 compound. 
In the previous studies KFe2Se 2 is thought to have the 
ThCr 2 Si 2 -type structure and K atom is located at the 2a 
Wyckoff position similar to Ba atom in BaFe2As2, shown 
in Fig. Qla). The 2a position is at the body center of 
the Se tetragonal formed by eight Se atoms between two 
neighbor Fe layers, while the 2b and 4c positions are at 
the middle point of the side edge and the center of side 
face in the Se tetragonal body, see Fig. Hfb) and (c). 

We perform the structure optimizations of KFe2Se2 in 
which the K atoms is sited in 2a, 2b and 4c Wyckoff posi- 
tion in space group I4/mmm respectively. The calculated 
results are showed in the Table. Uand the optimized lat- 
tice parameters c are compared with the measured ones. 



TABLE I: The lattice parameters c, the Fe-Se-Fe bond angle 
a and the Se height Hs e to Fe plane of KFe2Se2 in several 
magnetic orderings for the Struc-2a, Struc-2b and Struc-4c 
phases. In Bicoll-AFM order the lattice constants a shrinking 
lead to the decrease of a along the FM direction and the 
distortion of Fe layer lead to the fluctuation of Hse , shown in 
bracket. For the Struc-2b the Bicoll-AFM state is not stable. 



structure 


magnetic phase 


c(A) 


Hse(A) 


O 




Nonmag 


13.30 


1.29 


112.3 


Struc-2a 


FM 


14.04 


1.62 


98.7 




Neel-AFM 


13.67 


1.37 


109.5 




Coll-AFM 


13.79 


1.42 


107.8 




Bicoll-AFM 


14.45 


1.45 


107.0 








(1.64) 


(97.2) 


Exper. [lj 




14.04 








Nonmag 


18.60 


1.31 


110.6 


Struc-2b 


FM 


19.04 


1.59 


99.5 




Neel-AFM 


18.97 


1.39 


108.1 




Coll-AFM 


18.84 


1.40 


108.8 




Bicoll-AFM 








Exper. [8] 




20.48 








Nonmag 


15.46 


1.30 


113.9 


Struc-4c 


FM 


16.70 


1.57 


102.1 




Neel-AFM 


15.75 


1.35 


111.2 




Coll-AFM 


15.79 


1.38 


110.7 




Bicoll-AFM 


16.24 


1.42 


109.2 








(1.64) 


(96.7) 


Exper. [8] 




16.16 
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When K atom take up the 2a position, the relaxed lat- 
tice constant c is about 13.3 ~ 14.5 A for different mag- 
netic phases, consistent with the measured data 14.04 A 
synthesized by high-temperature routes. In the case of 
K occupying the 2b and 4c position, the optimized pa- 
rameters c are about 18.6 ~ 19.0 Aand 15.5 ~ 16.7 A 
respectively for different magnetic orderings, with good 
agreement with the observed lattice c 20.48 A and 16.16 
A synthesized using the ammonothermal method. 

The T c dependence of the As height to Fe plane and 
the Fe-As-Fe angle have been investigated extensively in 
iron arsenide superconductor and the neutron diffraction 
measurements reveal that T c becomes maximum when 
the Fe-As-Fe angle is close to 109.4 degree, correspond- 



ing to the perfect FeAs4 tetrahedron (16J. The calcu- 
lated Se heights to Fe plane Hs e and the Fe-As-Fe an- 
gles a in three structure phases of KFe2Se2 are presented 
in Table. HI In general, the Se heights within the given 
anti-ferromagnetic orders for three structure phases are 
in the range of 1.35 ~ 1.45 A and the the discrepancy 
of Fe-As-Fe angel between our calculation and the ideal 
value 109.4° is less than 2.4° (except the angles along 
ferromagnetic direction of Bicoll-AFM order in Struc-2a 
and Struc-4c). Since the Coll-AFM magnetic phase is 
the lowest or second lowest energy state for three struc- 
ture phases of KFe2Se2 (see Table ITT)) , it is convenient for 
us to compare the Fe-Se-Fe angles a in KFe2Se2 in this 
magnetic order. In our calculations the angles in Struc- 
2a, Struc-2b and Struc-4c phases are 107.0°, 108.8°, and 
110.7°. Among three structures, the 108.8° in Struc-2b 
is the most close to 109.4°. According to the qualitative 
conclusion that the ideal tetrahedral FeAs4 or FeSe4 cor- 
respond to the maximum of T c , we suggest that Struc-2b 
is related to the T c =40 K, while Struc-2a and Struc-4c 
correspond to T c about 30 K superconduting phases. As 
for the Se height to Fe plane, it is dependent on both the 
Fe-Se-Fe angle and the crystal lattice constants a (or 6), 
there is no definite correlation between Se height and T c . 

For the second part, we investigate the electronic prop- 
erties of the Struc-2a, Struc-2b and Struc-4c structures 
in nonmagnetic state. Fig. [H (a), (b) and (c) show the 
total and orbital resolved partial density of states (DOS) 
for three structures respectively. The DOS ranging from 
-7 to -3 eV are formed by the hybridization between Se 
4p and Fe 3d orbitals. The DOS at Fermi level is mainly 
from Fe 3d orbital and a small contribution is from Se 4p 
orbital. The Struc-2b has the maximum 3.37 states/eV 
per unit cell, compared with 2.15 of Struc-2a and 2.43 of 
Struc-2c, which reveal that Struc-2b should be possessed 
of some especial electronic properties different from the 
Struc-2a and Struc-4c phase. 

Figj3] (a), (d) and (c), (f) show the band structures 
and Fermi surfaces of the Struc-2a and Struc-4c phase 
respectively. The two structure phases have a great re- 
semblance in the band structure and the topological char- 
acter of Fermi surface, including that two Fermi surface 




-2 2 

Energy (eV) 

FIG. 2: (Color online) Total and orbital-resolved partial den- 
sity of states per primitive cell of KFe2Se2 for the Struc-2a, 
Struc-2b and Struc-4c structures. The Fermi level is set to 



sheets at the Brillouin Zone corner and a Fermi surface 
pocket around Z point. They also share the critical fea- 
ture that no Fermi surface appear near T point obviously 
different from the iron arsenide superconductor. Accord- 
ing to this resemblance, the Struc-2a and Struc-4c struc- 
ture should have the similar superconducting behaviors 
and the T c should be approximate to each other, associ- 
ated with the facts that T c are about 30 K in both the 
c=16.16 A and the c=14.04 A superconducting phases in 
KFe 2 Se2 compound. 

From FigO (b), (e) we can see that there are three 
bands crossing the Fermi level for the Struc-2b and the 
Fermi surfaces present better 2-dimensionality characters 
than that of Struc-2a phase. The most surprising thing 
is that there exist three cylinder shape Fermi surfaces 
centered around T-Z, the innermost is hole- type derived 
from the band 1 marked in Fig[3] (b) and the other two 
is electron- type derived from the energy band 2. The 
concentric cylinder Fermi surface sheets along T-Z is the 
typical features in iron arsenide superconductor, so the 
electronic structure of the Struc-2b phase of KFe2Se2 is 
very similar to the electronic structure of iron arsenide su- 
perconductor such as LaFeAsO and BaFe2 As2 . In Struc- 
2b phase, it is the emergence of cylinder Fermi surface 
around T-Z that is expected to bring about the stronger 
superconducting electron pairing than that in the Struc- 
2a and Struc-4c phases. The experimental T c = 40 K su- 
perconducting phase should correspond to the Struc-2b 
phase in KFe2Se2 compound. 

For the third part, we discuss the magnetic interaction 
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FIG. 3: (Color online) The energy band structures of KFe2Se2 
for the Struc-2a, Struc-2b and Struc-4c structure phase shown 
in (a), (b) and (c), and the corresponding Fermi surfaces 
shown in (d), (e) and (f) respectively. The bands crossing 
the Fermi energy denoted by line 1, 2 and 3 in (b). The 
Fermi level is set to zero. 



between Fe moments and its connection to the supercon- 
ducting temperature T c . The magnetic moments around 
each Fe atom were found to be about 2.8 [iB in the 
alkali metal intercalated iron selenide superconductors 
fH, much larger than the moments in the iron arsenide 
superconducting materials, so the magnetic interactions 
among Fe atoms play an important role in determination 
of structural and electronic properties of AFe2Se2 (A=K, 
Rb, Cs and Tl). We perform the calculation for the Struc- 
2a, Struc-2b and Struc-4c structures under the nonmag- 
netic order (NM), ferromagnetic order (FM), checker- 
board antiferromagnetic order (Neel-AFM) , collinear an- 
tiferromagnetic order (Coll-AFM), and bi-collinear anti- 
ferromagnetic order (Bicoll-AFM). The magnetic ground 
states for the Struc-2a and Struc-4c are the Bicoll-AFM 
order. The lowest energy state for the Struc-2b is the 
Coll-AFM order, since under the Bicoll-AFM order K 
atoms will migrate to 4c position and lead to the Struc- 
4c structure in our calculations. The energies of three 
structures under several magnetic orders relative to the 
nonmagnetic energy of the Struc-2a are listed in the Ta- 
ble, nn 

To qualify the magnetic interactions, we can assume 
that the energy differences among these magnetic orders 
mainly originate from the interactions between each two 
Fe moments with spin S. Then the frustrated Heisenberg 
model with the nearest , next nearest and third next 
nearest neighbor couplings Ji, J2 and J3 can be used to 
describe the interactions in KFe2Se2 system, namely, 



whereas {ij), <C ij 3> and (((ij))) denote the summa- 
tion over the nearest, next-nearest, and next-next-nearest 
neighbors, respectively. From the energy data in Table. 
ITT1 the superexchange coupling parameters Ji and J2 are 
shown in the rightmost column in Table. |TTf The detailed 
calculation is referred to Appendix of Ref. Il7l ) . 

By comparing these exchange coupling values, we no- 
tice that the Struc-2a and Struc-4c phases have the simi- 
lar features in the magnetic interactions, but the Struc-2b 
phase is somewhat different from them. Here we should 
pay close attention to the next neighbor superexchange 
interaction J2 in three structures. Theoretically, the two 
band model t — J\ — J2 indicates that the T,~ in iron based 
superconductor is proportional to the J2 [l8j and we also 
have drawn the same conclusion in the previous first prin- 
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From 



ciples studies on the iron based materials 
this point of view, we can obtain the results that Struc-2b 
phase has the higher superconducting transition temper- 
ature T c than Struc-2a and Struc-4c, and the latter two 
structures have the similar Tc values, because the next 
neighbor superexchange interaction J2 = 20.74 meV/S 2 
in Struc-2b is larger than J2 — 13.08 meV/S 2 in Struc-2a 
and J2 = 11.86 meV/S 2 in Struc-4c. The opinion that 
the higher T c superconduting phase occur in the Struc- 
2b structure is consistent with our analysis about both 
the Fermi surface features and the Fe-Se-Fe angle. The 
consistency is easy to understand because the superex- 
change interaction J2 between two next nearest neighbor 
Fe atoms is bridged by Se atom and has a direct connec- 



TABLE II: The energy per Fe atom, magnetic moment around 
Fe atom and the nearest and next nearest neighbor superex- 
change coupling of KFe2Se2 system in different magnetic or- 
derings for the Struc-2a, Struc-2b and Struc-4c structures. 
The nonmagnetic energy of Struc-2a is set to zero. For Struc- 
2b, the Bicoll-AFM state is not stable. 



H 





magnetic 


Energy 


moment 


Jl,J2 


structure 


phase 


(meV/Fe) 


( M B) 


(meV/S 2 ) 




Nonmag 


0.0 







Struc-2a 


FM 


-205.48 


3.20 


Ji=-25.27 




Neel-AFM 


-104.42 


2.25 


J 2 = 13.08 




Coll-AFM 


-207.41 


2.75 






Bicoll-AFM 


-248.25 


2.75 






Nonmag 


512.97 







Struc-2b 


FM 


388.02 


3.10 


Ji= 2.79 




Neel-AFM 


376.85 


2.36 


J 2 =20.75 




Coll-AFM 


298.32 


2.85 






Bicoll-AFM 










Nonmag 


307.83 


0.0 




Struc-4c 


FM 


77.39 


3.16 


Ji=-24.48 




Neel-AFM 


175.31 


2.37 


J 2 = 11.86 




Coll-AFM 


78.90 


2.79 






Bicoll-AFM 


30.77 


2.79 





«<y»> 
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tion to the Fe-Se-Fe angle. 

In addition, we do some discussion about the issue why 
K atom can locate in 2b or 4c Wyckoff position. Fixed 
the lattice parameters to the experimental values, we op- 
timize the K positions with initial K position sited in the 
2a, 26 or 4c respectively. The K positions do not change 
from the initial 2a or 4c to 2b position with the fixed c 
axis at 20.48 A, or from the initial 2a or 2b to 4c with the 
fixed c axis at 16.16 A. These calculations indicate that 
the Struc-2b and Struc-4c are the metastable structure 
phases relative to the Struc-2a. As to the formation of 
these metastable phases, we can assume the process hap- 
pened as blow. At first, the NH3 molecules and K atoms 
enter in between FeSe layers together and K atoms sit in 
2b or 4c Wyckoff positions, leading to the large length of 
lattice c. Then, NH3 escape from the compounds and K 
atoms still sit in 2b or 4c position. We also investigate the 
energy change with the lattice parameter c increasing for 
the Struc-2a, Struc-2b and Struc-4c respectively. When 
the length of lattice parameter c is large than 18.58 A, the 
energy of the Struc-4c begin to become less than the one 
of Struc-2a phase. When the parameter c is larger than 
20.85 A, the Struc-2b energy is the lowest one among 
three structure phases. The results mean that with the 
large lattice parameter c the K atom prefers to occupy 
the 4c or 2b Wyckoff position. 

In summary, we have investigated the effect of different 
Wyckoff position of K atom on the crystal structure of 
KFe2Se2 by the first principle calculations. We conclude 
that the occupations of Wyckoff position 2b and 4c by 
K atom correspond to the T c =40 K and T c =30 K su- 
perconducting phases in experiment respectively, which 
are demonstrated by the calculated lattice parameters, 
Fe-Se-Fe angle, electronic structures and superexchange 
coupling Ji- Our results have a good agreement with the 
experiment and provide a reasonable approach to under- 
stand the new superconducting phases in KFe2Se2- 
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